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Abstract 


A  general  framework  for  adaptive  and  non-adaptive  space-time  beampattem  synthesis  using  quadratic  beampattem 
constraints  with  linearly  constrained  minimum  variance  (LCMV)  beamforming  has  been  developed.  Main  beam 
and  sidelobe  pattern  control  is  achieved  by  imposing  a  set  of  inequality  constraints  on  the  weighted  mean-square 
error  between  the  adaptive  pattern  and  a  desired  beampattem  over  a  set  of  angle-Doppler  regions.  An  iterative 
procedure  for  satisfying  the  constraints  is  developed  which  can  be  applied  as  post-processing  to  standard  LCMV 
beamformers.  The  algorithm  is  used  to  synthesize  a  nearly  uniform  sidelobe  level  quiescent  pattern  for  the  circu¬ 
lar  UHF  Electronically  Scanned  Array  (UESA),  and  to  control  sidelobe  levels  for  the  same  array  in  an  adaptive 
manner.  The  technique  has  been  generalized  for  general  rank  reducing  transformations  to  reduce  computational 
complexity.  Performance  results  using  data  provided  by  Lincoln  Lab  show  that  under  low  sample  support  condi¬ 
tions,  sidelobes  can  be  effectively  suppressed  while  maintaining  high  signal-to-interference  plus  noise  ratio,  and 
deep  nulls  on  clutter  and  interferers. 
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1  Introduction 


Space-Time  Adaptive  Processing  (STAP)  used  in  airborne  radar  systems  combines  signals  from  N  antenna  array 
elements  and  M  pulses  to  adaptively  suppress  clutter  and  jamming  in  both  the  space  (angle)  and  time  (Doppler 
frequency)  dimension  [1].  Traditionally,  STAP  systems  have  used  a  rotating  linear  array  configuration,  however 
a  fixed  circular  ring  array  is  currently  under  development  under  the  UHF  Electronically  Scanned  Array  (UESA) 
program  sponsored  by  the  Office  of  Naval  Research  (ONR).  The  array  consists  of  54  directional  antenna  elements 
with  suppressed  backlobes.  Only  20  of  the  elements  will  be  used  at  a  time  to  transmit  and  receive  [2],  With 
this  configuration,  the  antenna  can  be  scanned  mechanically  in  6.67°  increments  by  choosing  the  appropriate 
20-element  sector,  and  scanned  electronically  ±3.33°  with  the  chosen  sector  of  elements.  The  circular  array 
configuration  has  the  potential  to  provide  continuous  360°  availability,  however  it  has  some  potentially  negative 
impacts  for  STAP  algorithms  [2],  First,  the  clutter  rank  is  increased  in  a  manner  similar  to  the  increase  from 
misalignment  with  the  velocity  vector  in  linear  arrays.  Second,  the  clutter  locus  varies  with  range.  This  decreases 
the  number  of  range  gates  that  can  be  averaged  to  reliably  estimate  the  clutter  covariance  matrix. 

The  foundation  of  most  STAP  techniques  is  the  Minimum  Variance  Distortionless  Response  (MVDR)  pro¬ 
cessor  [1].  The  standard  MVDR  processor  weights  are  designed  to  minimize  the  processor  output  power  subject 
to  a  linear  distortionless  constraint  in  the  angle-Doppler  steering  direction.  The  MVDR  beamformer  can  have 
unacceptably  large  sidelobes  and  mainlobe  squinting  due  to  sensor  perturbations,  pointing  error,  and  low  sample 
support.  In  radar  systems,  this  behavior  can  lead  to  increased  false  alarms  from  clutter  and  unexpected  interferes. 
Techniques  for  improved  pattern  control  include  steering  vector  tapering  using  a  quiescent  pattern  with  desirable 
main  beam  and  sidelobe  characteristics  [3],  linear  main  beam  constraints  [4],  white  noise  gain  constraints  [5],  and 
quadratic  quiescent  pattern  constraints  [6]-[10],  and  reduced  rank  subspace  techniques  [1 1]-[13]. 

To  mitigate  this  problem,  we  have  developed  a  general  framework  for  adaptive  and  non-adaptive  beampattem 
synthesis  for  non-linear  arrays  based  on  LCMV  beamforming  with  quadratic  beampattem  constraints  (QPC).  Some 
of  the  results  have  been  published  in  [14]-[17].  Additional  publications  are  in  preparation  (including  [18]-[19]). 
In  this  technique,  main  beam  and  sidelobe  pattern  control  is  achieved  by  imposing  a  set  of  inequality  constraints 
on  the  weighted  mean-square  error  between  the  adaptive  pattern  and  a  desired  beampattem  over  a  set  of  angle- 
Doppler  regions.  An  important  feature  of  the  LCMV-QPC  formulation  is  the  specification  of  multiple  quadratic 
pattern  constraints.  By  proper  choice  of  the  number  of  constraints,  the  angle-Doppler  regions  to  which  they  apply, 
and  the  desired  beampattems  in  those  regions,  the  level  of  pattern  control  can  be  traded  off  against  algorithmic 
complexity.  At  one  extreme,  low-complexity  techniques  can  be  obtained  based  on  one  or  two  constraints  similar 
to  the  adaptive  pattern  control  methods  in  [6]-[9].  At  the  other  extreme,  we  can  achieve  tight  pattern  control  using 
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many  constraints,  in  a  manner  similar  to  the  technique  in  [10].  Between  the  two  extremes,  our  approach  using 
several  constraints  has  been  shown  to  achieve  good  pattern  control  and  maintain  a  high  SINR  with  reasonable 
complexity.  The  algorithm  uses  an  iterative  procedure  for  satisfying  the  constraints  which  can  be  applied  as  post¬ 
processing  to  standard  STAP  processors.  The  technique  can  also  be  used  for  non-adaptive  pattern  synthesis.  It 
generalizes  the  techniques  in  [20]  and  [10]  for  developing  low  sidelobe  quiescent  patterns  for  arbitrary  arrays, 
and  can  be  used  for  developing  the  desired  patterns  used  in  the  adaptive  method.  In  this  report,  we  describe  the 
algorithm  and  present  circular  array  STAP  results  with  data  provided  by  MIT  Lincoln  Lab  [21]. 

2  LCMV  Beamforming  with  Quadratic  Pattern  Constraints 

2.1  Direct  LCMV  Problem  Formulation 

We  assume  a  STAP  model  with  N  antenna  elements  and  M  pulses.  Let  v(0,  p,  u>)  denote  NM  x  1  space-time 
array  response  vector  to  a  signal  arriving  with  elevation  angle  0,  azimuth  angle  p,  and  Doppler  frequency  ui. 
We  partition  azimuth  angle-Doppler  space  into  r  sectors,  fii, . . . ,  Or,  as  shown  in  Figure  1.  In  this  illustration, 
the  elevation  angle  space  has  only  one  partition  and  the  sectors  are  cubes,  however  more  general  partitions  of 
azimuth  angle,  elevation  angle,  and  Doppler  space  can  be  used.  Let  Bdti(0,  p,  u>)  =  w^v(0,  p,  u>)  be  a  desired 
beampattem  in  the  region  fij,  and  Wd,i  be  the  corresponding  weight  vector.  The  MSE  between  the  beampattem 
generated  by  the  adaptive  weight  vector  w  and  the  desired  beampattem  over  the  region  £1,  is  given  by 

e-  =  [  |wHv(0,  p,u>)  -  w^v^p,^!2^.  (1) 

J  fi  j 
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The  error  can  be  written  compactly  as 


ef  =  (w  -  wd,i)HQi(w  -  wd)i)  (2) 

where 

Qi  =  [  v(6,<f>,w)v(6,<t>,u))HdSl.  (3) 

Thus  the  pattern  error  is  a  quadratic  function  of  the  adaptive  weight  vector. 

Adaptive  weights  are  designed  according  to  the  standard  LCMV  criterion,  while  limiting  the  deviations  from 
the  desired  pattern  using  quadratic  pattern  constraints.  Let  R  be  the  data  covariance  matrix,  C  be  the  NM  x  d 
constraint  matrix,  and  f  be  the  d  x  1  vector  of  constraint  values.  The  LCMV-QPC  optimization  problem  is 

min  whRw  st.  CHw  =  f  (4) 

st.  (w  -  wd)i)HQi(w  -  wdii)  <  Li  i  =  l,...,r 

The  covariance  matrix  R  may  be  the  noise  and  interference  covariance  matrix  Rn  or  the  signal  plus  noise  and 
interference  covariance  matrix  Rx,  depending  on  the  application.  In  practice,  neither  one  is  known  exactly  and 
must  be  estimated  from  data.  In  many  radar  problems,  signal-free  training  data  is  used  to  estimate  Rn.  In  most 
other  applications,  including  active  sonar  and  communications,  the  signal  of  interest  is  always  present  and  Rx 
must  be  estimated.  One  commonly  used  estimate  is  the  Tv -sample  covariance  matrix 

Rx  =  i^x(fc)x(^  (5) 

where  x(k)  is  the  N  x  1  vector  of  array  data  at  time  k. 

The  LCMV-QPC  weight  vector  has  the  form 

w  =  Rq1C  ^C^Rq1^1  f  (6) 

+  Irq1  -  Rq'C  (c^Rq1^-1  ChRqX|  wq 

where 

r 

R  Q  =  R  +  ^AiQ,  (7) 

i—1 

r 

W  Q  =  (8) 

»=i 

This  is  the  multiple  constraint  extension  of  the  quadratically  constrained  MVDR  processor  developed  in  [8],  [9]. 
It  can  be  simplified  further  by  defining 

Pq  =  Rq1  —  RqXC  (c^R^c)-1  C^Rq1.  (9) 
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The  solution  then  becomes 


w  =  R^C  (c^R^c)  1f  +  PQwQ. 


(10) 


In  this  form,  we  see  that  the  first  term  satisfies  the  linear  constraint,  while  the  second  term  is  orthogonal  to  the 
constraint  vector  C  and  provides  additional  pattern  control. 

In  this  processor  a  weighted  sum  of  ‘loading’  matrices  Q,,  i  =  1 . . .  r  are  added  to  R,  and  a  weighted  sum 
of  desired  weight  vector  terms  QiV?d,i,i  =  1 . . .  r  appears  in  the  second  term.  The  loading  factors  balance  the 
adaptive  pattern  with  the  desired  pattern.  The  relative  contribution  of  these  terms  can  be  adjusted  to  achieve 
pattern  control  while  maintaining  high  signal-to-interference-plus-noise  ratio  (SINR).  There  are  generally  a  set  of 
optimum  loading  levels  A*,  i  =  1 . . .  r  that  satisfy  the  constraints,  however  there  is  no  closed  form  solution  for  the 
loading  levels,  even  when  r  =  1.  It  can  be  shown  that  the  mean-square  pattern  error  decreases  with  increasing 
but  at  the  expense  of  decreased  interference  suppression.  The  loading  levels  must  be  chosen  judiciously  to  achieve 
pattern  control  while  maintaining  high  signal-to-interference-plus-noise  ratio  (SINR). 

When  the  array  is  a  uniformly  spaced  linear  array,  and  there  is  one  sector  which  covers  the  entire  angle- 
Doppler  space,  Q  =  I,  and  the  technique  is  the  same  as  diagonal  loading.  When  one  or  two  quadratic  constraints 
are  imposed  over  the  main  beam  and/or  sidelobe  regions,  beamformers  similar  to  those  developed  in  [6]-[8]  are 
obtained.  When  there  is  a  single  main  beam  constraint  and  many  constraints  over  a  dense  grid  of  discrete  points 
in  the  sidelobe  region,  the  technique  in  [10]  is  obtained.  We  take  a  moderate  approach  using  several  constraints  to 
maintain  a  high  SINR  and  good  pattern  behavior  with  reasonable  complexity. 

2.2  Iterative  Implementation 

Here  we  present  an  iterative  update  procedure  for  the  partially  adaptive  LCMV-QPC  processor.  First  define 

q*  =  QjWrf,i.  (11) 

Now  the  quadratic  pattern  constraints  in  (5)  can  be  written  as 

wfiQ,w  -  2 5fte  (qfw)  <rji  i  =  l...r  (12) 

where 

r/i  =  Li  -  w^QjW^i.  (13) 

At  each  iteration,  the  pattern  errors  are  checked  against  the  constraints.  If  a  constraint  is  exceeded,  the  loading  for 
that  sector  is  increased  by  an  incremental  factor  A^p),  i.e.  A^p)  =  A^p_1)  +  A^p).  The  loaded  covariance  matrix 
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and  desired  weight  term  can  be  updated  according  to 


R?  =  R  +  E^Q* 


2=1 

r 

-  rJ-1)+Ea“c** 

2=1 

(14) 

=  I>Sp)q, 

2=1 

=  +EAiP)qi‘ 

2=1 

(15) 

The  adaptive  weights  can  then  be  computed  from 

S(p>  =  (Ref)1 

(16) 

p(p)  _  s (p)  -  S(p)c  (cHs(p)c)  1  CKS(p) 

(17) 

w(p)  =  s(p)c  (cHS(p)c)  1  f  +  P(p|w^. 

(18) 

This  procedure  isn’t  particularly  attractive  because  of  the  computational  complexity.  However,  if  the  incre¬ 
mental  loading  levels  A^p)  are  small,  the  update  can  be  accomplished  with  an  approximation  that  reduces  the 
complexity.  The  covariance  matrix  inverse  at  the  pth  iteration  can  be  expressed  as, 


s(p)  =  ns<p_1))  1  +  EAiP)Q 


-1 


(19) 


Substituting  (19)  and  (15)  into  (17)-(18),  we  can  expand  the  weight  vector  in  a  first  order  Taylor  series  approxi¬ 
mation  about  the  incremental  loading  levels  A \p\i  =  1 . . .  r.  After  some  fairly  tedious  manipulations,  this  gives 
a  direct  update  for  w(pl  and  P  ■p-1  as  follows: 


w(p)  _  w(p-i)  _  p(p-i)  E  Ajp)(QiW(p_1) 
2=1 


(20) 


p(p)  _  p(p-i)  _  p(p-i) 


EAiP)Q<  P(P_1)- 


(21) 


\i=l  / 

One  way  to  achieve  fast  convergence  while  ensuring  that  the  small  update  assumption  is  valid  is  to  let  A-p'1 
be  a  fraction  of  the  of  the  current  loading  value,  i.e.  A^p)  =  a\\p),  where  a  in  the  range  0.3  to  1  seems  to  work 
well.  This  requires  that  the  initial  loading  level  be  non-zero.  One  possibility  is  to  initialize  all  of  the  loading  levels 
to  some  small  value,  i.e.  Af )  =  A0,  i  =  1 . . .  r.  If  the  initial  loading  is  small  enough,  the  initial  weight  vector  is 
essentially  the  standard  adaptive  LCMV  weight  vector  given  in  (7). 
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To  summarize,  the  algorithm  is  initialized  by 


-i 


S(0) 

p(°) 


R  +  Ao^Q 

,  z=l 

s(°)  _  s(0)c  (cHS(0)c)_i  CHS(0) 
s(°)c  ^ +  P(0)  [a0^ 

\  i= 1 


Qi 


At  each  iteration,  the  weights  are  updated  by 


(22) 

(23) 

(24) 


1.  for*  =  l,...,r 

ifw(p-1)frQ.w(p-1)-  2^e  (qf  w^1))  >  rji 
then  =  qA|p_1'  ,  else  A^  =  0 
A(P)  =  A(P-1)  +  A(p) 

2. Q(p)=^AfQi  (25) 

i=l 

3.  q(p)  =  Q(p)w(p-1}  -  ^2  AiP)<li  (26) 

1=1 

4.  w(p)  =  W(P~1)  -  p(p-Dq(p)  (27) 

5_  p(p)  _  p(p-t)  _  p(p-t)Q(p)p(p-t)<  (28) 

The  LCMV-QPC  technique  can  be  used  for  non-adaptive  pattern  synthesis  by  letting 

R=  [  v(6,  </>,ai)v(0, 4>,w)HdCl,  (29) 


which  is  the  same  as  (3)  for  a  single  sector  covering  the  entire  angle-Doppler  space.  It  generalizes  the  techniques 
in  [10]  and  [20]  for  developing  low  sidelobe  quiescent  patterns  for  arbitrary  arrays,  and  can  be  used  for  developing 
a  tapered  steering  vector  for  use  in  the  adaptive  methods. 
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3  Examples 

In  the  MIT  Lincoln  Lab  data  set  [21],  there  are  IV  =  20  elements  and  M  —  18  pulses  with  a  300  Hz  pulse 
repetition  frequency.  The  UESA  spatial  beampattem  is  shown  in  Figure  2  and  the  space-time  beampattem  is 
shown  in  Figure  3. 

The  beam-pattern  is  steered  to  4>  =  0°  and  w  =  60  Hz  for  a  range  of  50  km,  which  corresponds  to  6  =  -10.5°. 
Angle-Doppler  space  was  partitioned  into  one  elevation  angle  sector  0  €  (—11°,  —2°),  11  azimuth  angle  sectors 
(j)  e  (-12°,  12°),  ±(12°,  30°),  ±(30°,  60°),  ±(60°,  100°),  ±(100°,  140°),  ±(140°,  180°),  and  5  Doppler  sectors 
oj  <=  (-30, 30),  ±(30, 90),  ±(90, 150)  Hz  for  a  total  of  1x11x5  =  55  sectors.  The  desired  pattern  was  set  to  zero 
outside  of  the  mainlobe  region,  and  the  constraint  levels  were  chosen  for  -35  dB  sidelobe  levels.  No  constraint 
was  used  in  the  mainlobe  region.  There  were  two  30  dB  interference-to-noise  ratio  (INR)  jammers  at  60°  and 
-20°,  in  addition  to  clutter.  An  8  km  training  window  (200  snapshots)  was  used  to  estimate  the  covariance  matrix. 
The  standard  LCMV  processor  weights  were  computed  by  adding  -30  dB  diagonal  loading  to  allow  the  covariance 
matrix  to  be  inverted.  The  resulting  space-time  beampattem,  and  beampattem  cuts  are  shown  in  Figures  4  and  5. 
The  beamformer  has  put  nulls  on  the  clutter  ridge  and  the  two  jammers,  however  the  sidelobes  are  quite  high. 

The  LCMV-QPC  adaptive  beamformer  was  used  to  reduce  the  sidelobes.  The  initial  loading  levels  were  set  to 
A0  =  0.0013,  and  then  iteratively  increased  using  a  =  0.8.  In  10  iterations,  the  LCMV-QPC  beamformer  is  able 
to  reduce  the  sidelobes  below  the  -35  dB  level  while  maintaining  a  well  behaved  main-beam,  and  deep  clutter  and 
jammer  nulls.  The  final  beampattem  is  shown  in  Figures  6  and  7. 
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Figure  5:  Initial  Adaptive  Beampattem 
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Iteration  10,  AG  =  17.9759  dB 
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Figure  6:  Final  Adaptive  Beampattem 


Beam  pattern  cut  at  60  Hz 


Figure  7:  Final  Adaptive  Beampattem 
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4  Refinements  and  Enhancements 

The  core  LCMV-QPC  technique  has  also  been  formulated  for  the  generalized  sidelobe  canceller  (GSC)  [22]  form 
of  the  LCMV  processor.  These  two  formulations  are  mathematically  equivalent  and  provide  flexibility  in  im¬ 
plementation.  They  are  based  on  an  earlier  technique  developed  under  ONR  Grant  #N00014-99- 1-0691  using  a 
minimum  mean  square  error  (MMSE)  approach  for  beamforming.  Although  the  LCMV  and  MMSE  approaches 
are  similar,  we  prefer  the  LCMV  approach  because  it  allows  specification  of  a  distortionless  constraint  in  the 
steering  direction  and  scaling  of  the  beampattem  is  not  an  issue. 

The  technique  has  been  generalized  for  general  rank  reducing  transformations.  This  is  important  because 
while  full-dimension  STAP  algorithms  provide  excellent  suppression  of  clutter  and  jammers,  they  require  a  large 
number  of  samples  («  2 NM)  for  optimum  performance  and  are  computationally  expensive.  Reduced-dimension, 
or  partially-adaptive,  techniques  operate  in  a  reduced-dimensional  subpace  by  exploiting  the  low-rank  proper¬ 
ties  of  the  clutter  and  interference.  Techniques  such  as  element-space  pre-Doppler,  element-space  post-Doppler, 
beamspace  pre-Doppler,  beamspace  post-Doppler,  PRI-staggered  post-Doppler,  and  eigenspace  have  been  shown 
to  reduce  computational  complexity  and  sample  support  requirements  while  providing  performance  close  to  full- 
dimension  STAP.  The  reduced  rank  LCMV-QPC  technique  can  be  applied  directly  to  these  techniques  to  provide 
additional  pattern  control. 
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